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We explore the effects of elastic scattering between dark matter and baryons on the 21-cm signal 
during the dark ages. In particular, we consider a dark-matter-baryon interaction with a cross 
section of the form a = aov~'^, in which case the effect of the drag force between the dark matter 
and baryon fluids grows with time. We show that, as opposed to what was previously thought, this 
effect heats up the baryons due to the relative velocity between dark matter and baryons. This 
creates an additional source of fluctuations, which can potentially make interactions easier to detect 
by 21-cm measurements than by using the cosmic microwave background and the Lyman-a forest. 

Our forecasts show that the magnitude of the cross section can be probed to (Tq ~ 3 x 10“'*^ cm^ 
for ^ 1 GeV and cro ~ 2 x 10”“^^ (m^^/lOGeV) cm^ for ^ 1 GeV with next generation 
experiments, and improved to (tq ~ 4 x cm^ for 1 GeV and cro ~ 4 x 10“'^^ GeV) 

cm^ for m-x 1 GeV with futuristic experiments. 


I. INTRODUCTION 

The standard picture of cold dark matter (CDM) [T] 
seems to fit very well with our current observational con¬ 
straints [2] . There are, however, a few puzzles that would 
require dark matter to have non-zero interactions [MB]. 
Moreover, several models for the dark-matter (DM) par¬ 
ticle predict some level of weak non-gravitational inter¬ 
action with standard-model baryons [Till]. Here we will 
study these interactions during the dark ages. 

The simplest way to observe these interactions would 
be through direct detection experiments, such as Dark- 
Side lallini, LUX m and XENONIOO [E]. These ex¬ 
periments are very sensitive to large dark-matter masses 
but cannot constrain interactions for DM masses below 
~ 10 GeV due to the small recoil of the nuclei in any 
given interaction. A different probe would consist of in¬ 
direct early-time effects of these interactions. One exam¬ 
ple would be the modification of the small-scale power 
spectrum, due to the drag induced in the DM by the in¬ 
teractions [13], which would be observable in the cosmic 
microwave background (CMB), as well as in Ly-a forest 
measurements. Another example is CMB spectral distor¬ 
tions, which would be created by the indirect coupling, 
through baryons, of dark matter and photons in the very 
early universe M- These last two probes require in¬ 
teractions to be relevant at early times, so they are not 
sensitive to all velocity dependences. Some models for 
dark-matter-baryon interactions may elude constraints 
because interactions get stronger at later times. We will 
focus on one of those models, in which the interaction 
cross section is parametrized by cr = aov~'^, one real¬ 
ization of which would be dark-matter milicharge [7]. To 
constrain interactions at later times, a useful probe is the 
21-cm line during the dark ages. 

The dark ages are the period following primordial re¬ 
combination and preceding the formation of the first lu¬ 
minous objects. During this cosmic era the only known 
observable is the redshifted Hydrogen hyperfine transi¬ 
tion, which traces the neutral Hydrogen density m- 
This observable has been proposed as a probe of non 


gaussianities OUT], as well as other effects that would 
modify the small-scale power spectrum |18j . Further¬ 
more, it has recently been proposed for the study of 
DM-baryon interactions m- We go beyond the analysis 
in Ref. [I9] by including the effect of relative velocities, 
which turns out to change the results significantly. 

Interactions between baryons and dark matter can be 
detected through their effect on the brightness temper¬ 
ature of the 21-cm line. This brightness temperature is 
proportional to the difference between the spin temper¬ 
ature of the neutral Hydrogen and the CMB tempera¬ 
ture. In the standard scenario the spin temperature is 
coupled to the baryon temperature during the redshift 
range z ^ 30 — 200. This creates a departure between 
spin temperature and CMB temperature. As shown in 
Ref. [TB], if the baryons are cooled down (by interacting 
with a colder fluid, like the dark matter) the spin tem¬ 
perature will be lower, modifying the overall brightness 
temperature. 

We emphasize that these interactions do not cause just 
cooling of the baryons, but also heating. In the usual pic¬ 
ture of interaction between two fluids, the warmer fluid 
will lose energy toward heating up the colder one, while 
there will be no energy transfer if both fluids have the 
same temperature. However, if there is a relative veloc¬ 
ity between the two fluids -dark matter and baryons in 
our case- there will be an additional friction term that 
will tend to damp this relative velocity. The kinetic en¬ 
ergy lost in this manner will induce heating in both fluids. 
The magnitude of this effect depends on the initial rela¬ 
tive velocity, which is given by a Gaussian variable with 
a (3D) variance of ^ 29 km/s at kinematic decoupling 
(z « 1010) [201121]. 

The brightness temperature will then acquire an ad¬ 
ditional spatial dependence, through the local variation 
of the relative velocities. Quantifying this effect, we find 
that during the dark ages it creates an additional con¬ 
tribution to the power spectrum of 21-cm temperature 
fluctuations, which can be more than an order of magni¬ 
tude bigger at large scales than the standard one, even 
for values of the cross section allowed by current CMB 
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studies [T3] . We study the detectability of this new signal 
with an SKA-inspired interferometeilj and with a more 
futuristic proposed experiment. We also study how the 
global signal changes due to interactions and discuss the 
prospects for experiments such as NenuFAIij^ 

This paper is organized as follows. In Section |ll] we 
derive the drag and heating terms and find their effects 
on the dark-matter and baryon temperatures. Later, in 
Section III we study how the change in the baryon tem¬ 
perature affects the signal of the 21-cm line during the 
dark ages. In Section we carry out a detectability 
analysis of this signal. We discuss some generalizations, 
as well as other possible effects of the interactions, in 
Section [V| before drawing our conclusions in Section [Vll 


(z = 1010) [H]. This rms value as well as the full power 
spectrum of 1 ^ 6,0 can be simply extracted from standard 
linear Boltzmann codes [Hum- 

Elastic interactions between fluids with a relative ve¬ 
locity will have two different effects. First, they will tend 
to decrease the relative velocity and achieve mechanical 
equilibrium, which in our scenario will manifest itself as 
a drag on the relative velocity m- Second, they will 
thermally couple the fluids, tending to equilibrate their 
temperatures. 

We start by calculating the drag on the relative veloc¬ 
ity. 


B. Drag term 


II. EVOLUTION OF INTERACTING DARK 
MATTER AND BARYON FLUIDS 

In this section we will study how the interactions be¬ 
tween DM and baryons change their temperatures. To 
do that we will have to calculate the drag on the rela¬ 
tive velocity due to interactions with baryons, as well as 
the heating effect on both fluids. Our results will rely on 
the current understanding of relative velocities, so let us 
start with a brief review. 


A. Velocities 


Throughout the text we consider cross sections 
parametrized as cr = aov~'^. First we analyze the velocity 
change due to the collision with a baryon with velocity 
Vf,. In the center-of-mass (CM) frame the initial velocity 
of the DM particle will be 

= (V, - V,) ( 2 ) 

mb + m^ 

and in an elastic collision the final velocity can be 
parametrized by the angle toward which it is scattered, 
so the final velocity of the dark matter particle is 

V^(CM)./ ^ ^ (3) 


In the standard cosmological evolution, dark matter 
starts collapsing as soon as matter-radiation equality is 
reached. Baryons, however, cannot cluster due to radia¬ 
tion pressure, until they decouple from the photon back¬ 
ground. This difference in their evolution history gen¬ 
erates a relative velocity between the two components. 
After the baryons and photons kinematically decouple, 
at redshift z « 1010, this velocity redshifts away, since 
the baryons experience infall into the DM gravitational 
wells. Ref. [ 20 ] first pointed out that relative velocities 
affect the formation of small-scale structure. Their effect 
on the standard power spectrum of 21-cm fluctuations in 
the dark ages was studied in Ref. [Hj. 

At kinematic decoupling, the relative velocities V^h = 
V^ — Vf, follow a Gaussian distribution, where V^ and 
Vfe are the DM and baryon bulk velocities. Then the dif¬ 
ferential probability of having an initial relative velocity 
V^f,_o is given by 

e- 3 V^i.,o/( 2 UD 
ViV ^ bfi ) = (-^^2 ) 3/2 ’ 

where the value of the (3D) width of this distribution 
is Vrms = 29 km/s ^ 10“"^ c at kinematic decoupling 


where h is a unit vector. The change in velocity in a 
single collision (which is Galilean invariant, and hence 
frame independent) is 


Av^ = 


mb 

mb -f m^^ 


|Vx-Vhl 



Vx - Vfa \ 
Wx-^b\J 


( 4 ) 


To calculate the full effect of the interactions we 
need to include the rate at which interactions happen, 
and average over the velocities of the fluid elements. 
The rate of interactions in a particular direction dh is 
da/dfi \ — -Vblnb, where (t(|vj^ — v^l) is the cross sec¬ 
tion as a function of the relative velocity, and Ub is the 
number density of baryons (targets). The time derivative 
of the DM bulk velocity will then be 

=nb J j d^Vbfb\v^-Vb\ J dn^Av^, (5) 


and we can perform the inner integral, by plugging 
Eq. Q into Eq. ([^ and realizing it has to be propor- 
tionalto the only direction (v^ — Vf,) inside the integral, 
to find 

[ d^Vbfbi^x - vt) |vx - Vi,| a, 
at mb + J J 

( 6 ) 

where we have defined the momentum-transfer cross sec¬ 
tion as 


da 

dcosO 


^ https://www.skatelescope.org/. 
2 http://nenufar.obs-nancay.fr. 


a 


d{cos6) 


(1 —COS 0 ). ( 7 ) 
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Alternatively, we could have calculated the drag on the 
baryon velocity, which is given by exchanging y o & in 
Eq. Q, so that dVb/dt = —{p^/pb)dV^/dt. The relative 
velomy between the two fluids will then evolve as 


dV 


Xb 


dt 


mi, + 


j d^Vxfx j 


d 'Ui,fi,{-v^ - vt) |v,^ - Vi,| cr, 

( 8 ) 

where we have defined Pm = Pb + Px- 

To calculate the two integrals over velocities we define 
two new variables and Vth, as 


Vm = 


rUv rrih 

—Vy H- -Vb 

Tx ^ Tb 


m 

J- X' 


X 


rub 

n 


and 


Vth = Vy - Vfc, 

so that the velocity distributions / factorize 


(9) 

( 10 ) 


d^Vxfx 


d^Vbfb = / d^vthfth / d^v^fm- ( 11 ) 


Nothing will depend on v^, so we can just integrate it 
out, leaving then only the integral of the relative veloc¬ 
ity Vth- The distribution function /th of this velocity is 
a Gaussian displaced from the origin by Vyj, and with 
thermal width given by the sums of the baryon and DM 
widths, Ty/my + Ti,/mb. The integral to calculate hence 
reduces to 

—^= -- f d^VthfthVthVthd{vth), ( 12 ) 

dt mb + m^ J 

Focusing on the case in which the interaction cross 
section is parametrized as (f = aov~'^, the drag term is 
given by 


D{Vxb) = 


dY-^b 

dt 


mb + my 


(13) 


where we have defined r = V^^b/uth, and — Tb/mb -I- 
T^lm-^, which is the variance of the thermal relative mo¬ 
tion of the two fluids. The function F{r) is determined 
as 


F{r) = erf 




(14) 


which grows with r from zero at r = 0 to one at r —>■ oo. 


C. Heating 

We now study the second effect that interactions have 
on the dark-matter and baryon fluids, namely heating. 
Interactions between two fluids (1 and 2) with different 
temperatures will tend to heat up the colder fluid (in our 
case the cold dark matter) at the expense of the energy 
of the warmer fluid, tending to equalize their temper¬ 
atures. The heating rate is usually proportional to the 


temperature difference (Ti — T 2 ). We will show here that, 
if there is a relative velocity between the two fluids, the 
heating rate will also include a friction term that will 
heat up both fluids, independently of their temperature 
difference. 

There is an intuitive reason to expect a heating term 
even for equal-temperature fluids, if two fluids with the 
same temperature collide with a relative velocity, and 
then equilibrate, this final relative velocity should van¬ 
ish. The kinetic energy would hence get transformed 
into a higher final temperature for both fluids, due to 
conservation of energy. 

Let us calculate the heating rate Qb of the baryons 
in their instantaneous rest frame, where the change in 
energy will directly give us the heat instead of having 
to add bulk motions. A baryon changes its energy in a 
collision by AEb = mt,vcM • Avb = —m^vcu ■ Ba, 
where vcm = (mb'Vb -k myVy)/(m{, + m^). The heating 
of the baryonic fluid per unit time is 


dQb 

dt 


rubPx 
(my -k mb) 


d^Vbfb 




X a (|vy - Vhl) |vy - Vb| [vcM • (Vfe 


Vx)] 


(15) 


where we have already integrated over outgoing angles 
dn using Eqs. Q and Q. 

We perform this integral in Appendix and find 


dQb ^ 2mbPxCroe '2 (Ty - Tb) ^ Py m^mb , 

dt (ruy -k mt,)^-\/^Mjj^ Pm my -k mb ^ ^ 


(16) 

The first term, in the r —0 limit, was derived in [13 BS], 
but here we also find the second term, which is non-zero 
for r 7 k 0 . 

By symmetry, Qy is obtained by simply substituting 
& O X in Eq. (16). We see that these expressions, with 
the drag D(Ey^in Eq. (13), conserve the total kinetic 
energy density in the baryon-DM fluid, i.e. 



Ub- 


dQb 

dt 


^D{V^b)Vxb = 0 . 


(17) 


Now that we know how the interactions change the 
energy of the baryons and DM at any given time, let us 
find how their temperatures are modified. 


D. Temperature evolution 


Using the expressions for the drag D(V^b), in Eq. (13), 
and the heating rates Qb and Qy, in Eq. (16), we can 
write the equations of the temperature evolution, miis]. 
In our analysis we also evolve the relative velocity Uyf,. 
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The set of equations we will have to solve is then 


da 

n'd'x 1 

a 3aH' 

(18) 

d/Tb _ 
da 

^Tb be . . 2Qb 

(19) 

dVxb 

D{Vxb) 

( 20 ) 

da 

a aH ’ 


where we have assumed the photon temperature is 
unaltered, H is the Hubble parameter and Tc is the 
Compton interaction rate, which depends on the free- 
electron density Ue- Since the free-electron abundance 
also depends on the baryon temperature through the re¬ 
combination rate, we must solve for Eqs. (18)-(20) simul¬ 
taneously with the free-electron fraction = nejnH 


( 21 ) 

where C is the Peebles factor [53], Eq is the ground en¬ 
ergy of Hydrogen, and AsiTb^T^) and BsiT-y) are the 
effective recombination coefficient and the effective pho¬ 
toionization rate to and from the excited state respec¬ 
tively 1531155] , 

For convenience, we parametrize the results in terms 
of a dimensionless cross section 1741 , defined as 


= ^0 


( 22 ) 


so that tT 4 i < 3.2(rn;^/GeV) is the 95% C.L. constraint 
from CMB-analysis [13], valid only for ^ mt- 


E. Limiting cases 


To gain understanding of the implications of Eq. (16 1 it 


is enlightening to study the extreme cases of very-heavy 
and very-light dark matter. 


For very massive dark matter (m^ ^ uib 


1 


GeV), the first term in Eq. (16 1 is small and the sec¬ 


ond one dominates, which means that the new effect 
we have calculated is more relevant than the previously- 
known result. In this limit we then have Qb = 
{Px/Pm)'mbVyf-bD{V^b)[^+0{mb/m^)], which means Qb oc 
aolm^. Equivalently, the DM heating term will be given 
by Qx = {Pb/PTn)mbVxbD{Vxb)[^ + 0{mb/mx)], so that 
Qx oc aojmx as well, so for ^ mb the constraints we 
will find will behave as (Tq oc 

• In the opposite limit, in which <C mb, we find 
that the temperature-independent heating term (second 
term in Eq. is linear in m^ and hence subdomi¬ 

nant. The first term is roughly constant. Although itth 
depends on T^/m^, starts as zero and does not change 
unless there are interactions. This leads to a net mass- 
independent cooling Qb < 0 , whereas the dark matter 
decouples, since Qx oc m^^ —>■ 0 . 

Let us now briefly discuss the two limiting cases where 
either thermal or relative velocities dominate. 


• When Vxb < Wth = V%7 mx + Tb/mb (thermal ve¬ 
locity dominates), we recover the results of Ref. [19], 
where baryons get cooled down and tend to thermalize 
with the dark matter fluid. This is shown in Fig. [^ as 
the “ 14 ;b,o = 0” case. 


• In the limit where Vxb is much bigger than uth, the 
second term in Eq. (16) dominates, which causes a net 
heating of the baryon fluid. However, the overall rate of 
interactions (and hence net heating or cooling) is sup¬ 
pressed for large velocities, due to the fact that the cross 
section is proportional to v~'^. 


F. Numerical results 


We solve the system Eqs. (18)-(21) for different values 


of (T41 and mx, starting at z = 1010 with the baryons 
tightly coupled to the photon fluid {Tb — T^) and with 
perfectly cold dark matter {T^ = 0 ), although we tested 
that having slightly warm dark matter at recombination 
does not change our results significantly. We use cosmo¬ 
logical parameters consistent with their current best-fit 
values [5]. We have also checked that, for the values of 
CT 41 considered in our analysis, the system is not already 
tightly coupled at z = 1010 , which would require us to 
start evolving the system at an earlier redshift. 

As for the initial conditions for Vxb, we will solve the 
system for an array of values from zero initial velocity to 
three times the width of its Gaussian distribution. For 
purposes of illustration we will plot two different cases, 
one in which Vxb,o = Arms = 29 km/s at initial redshift, 
and another in which V^b,o = Oj to show how the relative 
velocity affects the results. In the case with Vxb,o ^ 0, 
higher values of imply a more significant heating of 
the baryons. 

In Fig.[l]we show how the baryon temperature changes 
with the strength of the interactions. In the central and 
bottom panels we have mx > vn^. In those two figures 
it is explicit that having Vxb,o ^ 0 (red lines) induces 
extra heat in the system as a result of the damping of the 
relative velocity, which increases the temperature of both 
baryons and dark matter. However, when considering 
the case with Vxb,o = 0 (blue lines), the interactions cool 
down the baryons and only heat up the dark matter. In 
the upper panel of Fig. we have set mx = 0.1 GeV. In 
this case it is clear that introducing interactions can only 
cool down the baryons, albeit with a more pronounced 
temperature drop in the V^b,o = 0 case. 


III. EFFECTS ON THE 21-CM DARK AGES 
SIGNAL 

We have seen how the baryon and dark matter tem¬ 
peratures change when adding interactions. Now we will 
study how this modified baryon temperature gives rise 
to a different spin temperature for the gas during the 
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1+z 



1+z 



1+z 


FIG. 1. Baryon temperatures (three upper curves) without 
interactions (solid curve) and when adding interactions with 
(T41 = 1 (dashed-blue curve for the case where = 0 and 

red curve for 14 . 6,0 = Irms), as well as dark-matter tempera¬ 
tures (two lower curves, dash—dotted-blue curve for the case 
where I4i,,o = 0 and red curve for = Kras). From top to 
bottom we show the results for = 0.1, 1, and 10 GeV. 


dark ages, which in turn modifies the 21 -cm brightness 
temperature we would observe. 


A. 21-cm brightness temperature 

The electronic ground state of neutral Hydrogen is split 
into two hyperfine states, a singlet spin -0 and a triplet 
spin-1 configuration. The singlet state has a smaller en¬ 


ergy, with the transition from the triplet to the singlet 
corresponding to a wavelength of 21 cm. Because of its 
very long wavelength it is hard to confuse with any other 
redshifted line, making it a very unique probe of the 
physics of the early universe m- 

We define the spin temperature of the baryon gas 
through the ratio of the populations of the triplet to the 
singlet states, 


’!b = 3e-^*/^4 (23) 

no 

where T* = 0.068 A = 5.9/reV is the energy correspond¬ 
ing to the 21-cm transition. During the dark ages the 
upward and downward transitions are much faster than 
the evolution of the universe. This means that we can 
use the quasi-steady-state approximation and, to a good 
accuracy, find the values of ni and uq for which there is 
equilibrium, 

no {Coi + Roi) = ni {C\o + i?io); (24) 

where Rij are the rates of radiative transitions of the 
CMB blackbody photons and Cij are the collisional tran¬ 
sition rates [15]. We will always have T* ^ in 

which case the spin temperature is very well approxi¬ 
mated by 


Ts=T^ + 


Cio 

Cio + ^ 10 ^ 


(25) 


where Aio is the downward spontaneous Einstein co¬ 
efficient of the 21-cm transition. We neglect the 
Wouthuysen-Field effect [2M28] that would arise from 
inelastic scattering of Lyman-a photons after the first 
stars are created. 

The fact that the spin temperature is smaller than the 
CMB temperature during the dark ages enables the Hy¬ 
drogen gas to resonantly absorb more CMB photons with 
a rest wavelength of 21 cm than it emits. We would ob¬ 
serve this effect as a smaller brightness temperature of 
the CMB. Let us define the temperature (T 21 ) of the 21- 
cm line (redshifted to today) as the brightness temper¬ 
ature contrast with respect to the CMB, which is given 
by 


T21 


Ts-T^ 
1 + z 



Ts-T^ 
1 + z 


(26) 


where we have used the fact that the optical depth r is 
very small, and is in fact given by 


r = 


3 T* ,3 Aio 

H{z) + {l + z)drVr' 


(27) 


where H{z) is the Hubble rate, uhi is the density of neu¬ 
tral Hydrogen, A 21 ~ 21 cm and drVr is the proper gra¬ 
dient of the peculiar velocity along the radial direction. 

In the usual scenario the spin temperature follows the 
gas temperature as of decoupling and until z ~ 30, which 
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makes it different from the CMB temperature in the red- 
shift range z ^ 30 to 200. This creates a non-zero 21-cm 
line temperature T 21 in this range. As we have shown, 
dark-matter-baryon interactions can either cool down or 
heat up the baryons, thus changing the spin temperature. 


We show this effect in Fig. where we plot for ref¬ 
erence the CMB temperature, as well as the usual non¬ 
interacting gas and spin temperatures. We also plot the 
gas and spin temperature for interacting cases with ei¬ 
ther V^b,o = 0 or V^b.o = Crms- The deviation of the 
spin temperature in the interacting cases (blue and red 
curves) is apparent, even for a cross section of (T41 = 1, 
compatible with CMB bounds. 

If there is more heating than cooling of the baryons, the 
21 -cm brightness temperature decreases in magnitude, 
since the spin temperature is closer to the CMB tem¬ 
perature during the dark ages. Cooling of the baryons 
increases the brightness temperature, as long as the spin 
temperature stays coupled to the baryons. In Fig. [3| w e 
plot the 21-cm brightness temperature, from Eq. (261, 
for different values of the relative velocity and DM mass. 
It is interesting to note that the heating increases with 
the mass of the dark matter, as predicted, so that the 
average brightness temperature T 21 during the dark ages 
is higher when including interactions. 


B. Global signal 

Let us define T 2 i{V^b) as the brightness temperature 
in the absence of density perturbations. In the standard 
scenario this quantity is spatially homogeneous and is 
termed the global 21-cm signal. Once DM-baryon inter¬ 
actions are included, T 2 i{Vy^b) is still a function of the 
initial relative velocities. We calculate its average over 
said initial velocities as 

(T 21 ) = J d^V^b,of2iiV^b,o)ViV^b,o), (28) 


with the probability distribution V^V^bo) given by 
Eq. (§. We show this quantity in Fig. for the interact¬ 
ing case and for three different DM masses. 


C. 21-cm fluctuations 



1 +z 



1 +z 



1 +z 


FIG. 2. Values of the spin temperature (dashed curves) and 
the gas temperature (solid curves) for the collisionless case 
(black curve) and when including collisions (blue curve for 
Vxb,o = 0 and red curve for V^bfi = E-ms), as well as the CMB 
temperature in the dashed-green curve. From top to bottom 
we show the results for = 0.1, 1, and 10 GeV. 


As we have shown, the brightness temperature T 21 of 
the 21-cm line is modified by the inclusion of interactions, 
and this modification depends on the initial relative ve¬ 
locity. The large-scale fluctuations of the relative veloc¬ 
ity will therefore be imprinted on the brightness tem¬ 
perature, since two regions with different initial relative 
velocities will appear with different brightness tempera¬ 
tures (compare blue and red lines in Fig. [^, which will 
actually generate an additional contribution to the power 
spectrum of the 21-cm fluctuations. Let us calculate it. 


The standard deviation of T 21 as a function of V^b,o is 

T21,rms = ^(T|i)-(r 2 i)". (29) 

Even if T 21 had no explicit spatial dependence, it would 
fluctuate because relative velocities are not homogeneous. 
In principle, to compute the power spectrum of T 21 , one 
should first compute its two-point correlation function. 
This is obtained by integrating over the six-dimensional 
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FIG. 3. Values of the average brightness temperature of 
the 21 -cm line for the collisionless case (solid-black curve), 
the case with interactions (blue-dashed curve for V^i,,o = 0 , 
purple-dot-dashed curve for = Kms), and the average 

over initial velocities in the red-dotted curve. From top to 
bottom we show the results for = 0.1, 1, and 10 GeV. 


joint probability distribution of the relative velocities at 
two different points (see Ref. [21] )• To simplify matters 
we shall make the following approximation 



1+z 

FIG. 4. Values of the brightness temperature of the 21-cm 
line for the collisionless case in solid-black curve, and three 
with interactions ( 0-41 = 1 ), in dashed-blue curve = 0.1 
GeV, in dot-dashed—purple = 1 GeV, and in dotted-red 

curve = 10 GeV. 


variance of T 21 . For illustration we show T 21 as a function 
of T^b.o for the = 1 GeV case in Fig. We calculate 
the power spectrum of T 2 i(K^b^o) in this approximation 
to be 


(r2i(k)r*i(k')) = Tl^,^,Py2^{k){27T)HD{k + k'), (31) 

where Py 2 is the power spectrum of •y/2/3(l — 
Vxb,o/Kl,s)- We plot Pv^^ik) in Fig. 
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FIG. 5. Brightness temperature T 21 of the 21-cm line for 
= 1 GeV and (J 41 = 1 at redshifts 2 = 30 in dashed- 
blue curve, 2 = 40 in dash-dotted-purple curve and 2 = 50 
in dotted-red. We also show the average over velocities for 
each redshift, as defined in Eq. (28 1 , in solid curves and their 
corresponding colors. 


T2i{V^b,o) 


(T 21 ) + T21,rms 




(30) 


which has the advantage of resulting in simple analytic 
expressions (22) while still reproducing adequately the 


Our observable, the brightness temperature of the 21- 
cm line, varies in space through its dependence on the 
baryon density Ub, as well as on the initial relative veloc¬ 
ities V^b,o- To linear order in density perturbations the 
temperature of the 21-cm line, Eq. (26), will be given by 







































Power Spectrum 


m 


T21=f2l{V^b,o) + ^S, 
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where k = /r"^ + A:j|, and W{k\\) is the window func¬ 
tion. The new angular power spectrum due to 

(32) interactions, will be 


where S = (rib —fib)/fib, dT 2 i/d 5 is a well-known function 
of redshift for V^b,o = 0, and T 21 = f(Ts — T^)/(l + z) 
is the unperturbed value of the brightness temperature. 
Both T 21 and dT 2 i/d 5 depend on the initial relative ve¬ 
locities. The average over initial velocities of T 21 is then 

{T 2 i) = {f 2 i) + (^^)5. (33) 

We can, however, approximate (dT 2 i/d 5 ) « 

dT 2 i/d 6 {V^b,o = 0)j since the error made in the 21-cm 
temperature would be of order 6 T 2 i(k) ^ T 2 i^rms^, which 
is subdominant. We can calculate the variance of T 21 
over both initial relative velocities and overdensities (as 
in the usual power spectrum) to find 


Pt 2 i (k) = (k) + f a(z) -I- 121 ^ j Pb(k, z), 

(34) 

where fc|| is the magnitude of k in the line-of-sight di¬ 
rection, a(z) as defined in Ref. and Pb is the usual 
baryon power spectrum. 



FIG. 6 . Power spectra as a function of k. In solid-blue c urv e 
we show the power spectrum of V^i,, as defined in Eq. (311, 
and in dashed-black curve the standard power spectrum of 
baryon overdensities at redshift 2 ; = 30. 


We can convert easily from A:-space to £-space by using 
a harmonic transform |3()j . which is exact in the case of 
the flat-sky limit and still a very good approximation 
for i > 10, which should be good enough for our order- 
of-magnitude estimates. We define the angular power 
spectrum for the standard fluctuations (std) as 


C std _ 

e — 


1 

j,2 


dk\\ 

2ti 


2 


W(k\\) 


/ n2 


Pb(k), (35) 


rp2 

^Vxb _ ^21,rnis 

^ ^12 


dfc|| 

27r 


W(fc||) 


P 


y 2 ^ (k) 


(36) 


Before going into a full-scale analysis one might be 
interested in what would happen at a single £, and at 
different redshifts. We show in Fig. the value of the 
square root of the velocity power spectrum 
for different values of the cross section. For illustra¬ 
tion purposes we also show the usual power spectrum 
(Cf‘'^)^P, from Eq. (35), where for simplicity we have 
taken a{z)/ a(zo)(Cf^P ( zq) as a proxy for {Cf‘^)^P(z) 
as well as a redshift-independent bandwidth of ^vjv = 
0.02 to avoid recalculating the integral in Eq. ( [^ for 
each redshift in this plot. We show the cases of ^ = 30 
and £ = 1000 in the upper and lower panels, respectively. 



1+z 



1+z 


FIG. 7. Amplitude of the brightness-temperature fluctuations 
induced by relative velocity fluctuations for three different 
cross sections (solid-black curve for (T 41 = 0 . 01 , red curve for 
CT 41 = 0 . 1 , and green curve for (T 41 = 1 ), as well as for the usual 
baryon density perturbations. We calculate at two different 
scales: in the first panel we show the results for f = 30 and 
in the second one for i = 1000 . 
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IV. DETECTABILITY 


So far we have shown that DM-baryons interactions 
modify the baryon temperature, raising it or lowering 
it, depending on the initial relative velocity. Varying 
the baryon temperature will change the spin temperature 
and hence the brightness temperature of the 21 -cm line. 
This quantity, also known as the “global signal”, is the 
main observable during the dark ages. We will study how 
to detect interactions with a global-signal experiment. 

Moreover, since the temperatures depend on initial ve¬ 
locities, and these have a spatial-dependence, we have ar¬ 
gued that there will be a new contribution to the power 
spectrum, which, at large scales, can overcome the stan¬ 
dard one for values of the cross section of 0-41 > 0.1. We 
will study the detectability of this signal with interfer¬ 
ometry later in this section. 


A. Global signal 


Let us start by analyzing the most direct effect of DM- 
baryon interactions, the change in the global signal dur¬ 
ing the dark ages. Next-generation experiments, such as 
NenuFAR, will survey the 21-cm line brightness temper¬ 
ature down to frequencies possibly as low as ^ 10 MHz, 
which corresponds to a redshift z > 100 . 

We have seen in Fig. how the brightness tempera¬ 
ture changes when adding interactions. We will use the 
amplitude of the brightness temperature at its peak as a 
proxy for the detectability of the signal, even though its 
very high redshift (z ~ 90) may make it unobservable. 

Let us first find the signal-to-noise ratio to detect inter¬ 
actions having a cross section CT 41 = 1. If we could deter¬ 
mine the brightness temperature T 21 at its peak with 5% 
precision, we would be able to detect interactions with 
tT 4 i = 1 at a signal-to-noise ratio S/N ~ 10 for = 0.1 
GeV, S/N ^ 0.5 for rriy = 1 GeV, and S/N > 1 for 
my = 10 GeV. 

More interestingly, if we were able to improve the 
error by a factor of 5, reaching 1% precision of peak- 
temperature determination, we would be able to detect 
cross sections as small as 0-41 < 0.04 for =0.1 GeV, 
0-41 < 0.1 for my = 1 GeV, and 0-41 < 0.2 for my = 10 
GeV, all of which are beyond what can be achieved by 
current GMB analysis m- 


B. Fluctuations 


We now turn our focus to the measurement of the 21- 
cm power spectrum, Eq. (34). In a maximum-likelihood 
analysis, the Fisher forecast for the error in the measure¬ 
ment of the amplitude A of a power spectrum Ci is given 


by [310 



(37) 


For a given sky coverage /sky, the error for an individual 
£ in the estimated value A is 1311131135] 


af = 


/sky(2^+l) 


{C\ 


'std |_ 


(38) 


where Cf is the ins tru mental noise power spectrum, de¬ 
fined below in Eq. (40), and is the standard power 
spectrum of 21 -cm fluctuations (under the null hypothe¬ 
sis of no DM-baryon interactions), Eq. (34). 

The minimum detectable amplitude A at l-tr signifi¬ 
cance is thus 


a 


A 


/sky ^ 

2 ^ 

^min 


(2£^1)(G^1 ' 

{c/ + c^Y _ 


(39) 


where = Cj^'’/A encodes the £ dependence of the 
velocity power spectrum from DM-baryon interactions, 
and £min = 180/0 is the largest scale accessible by an 
experiment with sky coverage /sky = 0^. Because of the 
use of the harmonic transform in Eq. ( |3^ , we take £yain = 
15. This should not affect the results significantly since 
there are very few modes at lower £. 

We will consider two different scenarios, first a realistic 
experiment modeled after SKA that could be taking data 
within the next few years and second a more idealized ex¬ 
periment whose noise level will be low enough to detect 
the primordial power spectrum at redshift z = 30 (but 
still not cosmic-variance limited, since the usual primor¬ 
dial power spectrum vanishes for smaller redshifts but 
the noise will not). 

We will study the redshift range z = 20 to 30, at the 
very end of the dark ages. This range is chosen to avoid 
complex astrophysical processes at low redshift as well as 
to still be observable from Earth. We may not be fully 
free of contamination, however, since the epoch of the 
formation of the first stars is unknown, and the X-rays 
generated during star formation may start to heat up 
the gas at z < 25 [36]. Moreover, accreting intermediate- 
mass black holes (sometimes termed miniquasars) may 
also be an important source of X-rays during this era 
[37H33] . Once data of the gas temperature during the 
dark ages are acquired, a careful analysis should take 
these processes into account along with the heating pro¬ 
duced by DM-baryon interactions, and by studying their 
different redshift behaviors and angular structures, dis¬ 
entangle them. We motivate future work to address this 
issue. 


® We assume that the likelihood function is Gaussian in the vicinity 
of its maximum 13211331 . 
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The angular noise power spectrum of an interferometer 
is given by [40] 


^^^ofcover \^cover(^)J 


(40) 


where ^cover(j^) = 2 ttDI\ is the maximum multipole at 
frequency v (corresponding to wavelength A) that can be 
measured with an array of dishes with maximum baseline 
D, covering a total area ^totai with a covering fraction 
/cover — ^dish^dishMtotai, in a frequency window A:/ 
with an observing time tg. The system temperature is 
given by Tgys 180 (i^/180 MHz)~ ^ K, consistent with 

un¬ 
inspired by design plans for the Square Kilometer Ar¬ 
ray, we first consider a future ground-based interferom¬ 
eter with access to the final stages of the dark ages, 
z CO 20 — 30, with a baseline of Z? = 6 km [corresponding 
to a maximum angular scale •Zcover('-^) ~ 5800 at redshift 
z = 30], with /cover = 0.02, surveying a sky fraction 
/sky = 0.75 for a total of five whole years. As for the 
bandwidth, we surveyed a range between Az/ = 0.1 MHz 
and 10 MHz and found that Av ~ 1 MHz is the opti¬ 
mum value (for smaller bandwidths the noise CiS dom¬ 
inate over the signal and for larger ones the number of 
redshift slices is too small). 

For more optimistic constraints, we set D = 50 km, 
/cover = 0 . 1 , and assume ten whole years of observations. 
In order to get a result closer to the cosmic-variance limit 
we could perform the analysis from z = 20 , going up to 
the beginning of the dark ages, z = 200. However, we find 
that it does not improve the results signihcantly, due to 
the rise of synchrotron radiation at low frequencies, which 
grows much more rapidly than the signal. We consider 
then the same redshift range as before, z from 20 to 30. 

One of the great advantages of 21-cm as a probe is 
the ability to analyze the tomography of the signal, en¬ 
abling us to coadd information from different redshift 
slices. Summing over redshift slices, the signal-to-noise 
ratio is given by 


(S/N) = 


(cy-k/ 


(41) 


In Fig. [^we show the C^s for the usual primordial per¬ 
turbations (Cl**^), for the instrumental noise {C^, both 
with next-generation and futuristic parameters), and for 
the new contribution due to interactions all of 

them at redshift z = 30. 


Results 

Let us start by considering the realistic noise case (that 
corresponds to the experimental parameters of SKA) and 
find what the signal-to-noise ratio would be for detect¬ 
ing (741 = 1. We calculate the signal-to-noise ratio for 



FIG. 8 . Angular power spectra at redshift z = 30 with band¬ 
width Au — 1 MHz. In solid-black curve we show the usual 
primordial perturbations, in solid- and dotted-blue curves the 
instrumental noises for the realistic and optimistic cases [see 
Eq. (401 and discussion below] and in dashed-red curve the 
new piece due to interactions for (J 41 = 1 and rriy = 1 GeV. 


CT 41 = 1 in each redshift bin between z = 20 and z = 30 
with Eq. ( jdTj ). We find the total signal-to-noise ratio to 
be S/N~ 3 for the case of rriy^ = 0.1 GeV, S/N~ 9 for 
rriy^ = 1 GeV, and S/N^ 0.2 for = 10 GeV. We could 
alternatively express the results in terms of the smallest 
CT 41 that would still give us a signal-to- noise ratio of 1 , 
taken to be approximately 1741 ,min = 1/-^ S/N. We show 
the minimum detectable cross sections in Tab. [ij 

Let us now move on to trying to find the smallest pos¬ 
sible 0-41 detectable at S'/ZV = 1 in the more optimistic 
case. In principle the amplitude A of , equal to 
T|i is a non-trivial function of redshift and < 741 . How¬ 
ever, we find that for small values of 1741 ( 0-41 < 0 . 1 ), the 
quantity /(z) = T 2 i_rms/i 74 i is approximately indepen¬ 
dent of cr 4 i (although it does depend on rriy). Then we 
can construct an estimator for an for each redshift slice. 



with variance given by ) = a\{z)jf'^{z). We can 
then combine all the estimators into a minimum-variance 
one, finding the variance of the final redshift-independent 
estimator. 





(43) 


With the optimistic experimental parameters defined 
above we find that the minimum cr 4 i observable at 68 % 
C.L. (Icr) is tT 4 i < 1.7 X 10“^ for = 0.1 GeV, 0-41 < 
4.3x10“^ for = 1 GeV, and 1741 < 3.6x10“^ for = 
10 GeV. These results are about 2 orders of magnitude 
better than the CMB constraints found in na, where 
cr 4 i ^ 16(mj^/10GeV). 
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nix [GeV] 

1/10 

1 

10 

Fluctuations (realistic) 
Fluctuations (optimistic) 
Global signal (1% error) 

6x10“*^ 

2x10“** 

4x10“*® 

3x10“*® 

4 X 10“** 
1x10“*® 

2x10“** 

4 X 10“*® 
2 X 10“*® 


TABLE I. Minimum ag (in cm^, corresponding to 0-41 x lO'^^) 
detectable with both realistic and optimistic interferometer 
parameters at 68 % C.L., as well as with global-signal analysis 
with 1 % accuracy for three different dark-matter masses 
(in GeV). 


V. DISCUSSION 



Before concluding we would like to make a few remarks: 

• As we have shown, interactions between dark matter 
and baryons give rise to a new heating term, which can 
increase the temperature of the baryons significantly. We 
only used that heating to study dark-ages physics but 
this result may have applications beyond our analysis, 
for example in the epoch of reionization |42l |43] . 

• In this work we have focused only on the case where 
(j v'^ with n = —4, but one may wonder whether 
the dark ages can potentially provide new information 
not contained in the CMB analysis for other values of n. 
Since the dark ages occur more recently than decoupling, 
we have only been interested in interactions that increase 
at later times. Ref. m showed that the interaction rate 
grows for n < —3, so all results that we could forecast 
for n > —3 would be worse than those obtained with 
CMB studies. That still leaves n = —3 as a potential 
interaction to study, for example. 

• It is also worth mentioning that if we wanted to 
translate these results to a constraint specific to a dark- 
matter milicharge model [7] , the ionization fraction of the 
baryons would cause a suppression of Xe ^ 10“'*. 

• We have also found a decrease in the bulk relative 
velocity of baryons and dark matter characterized by a 
drag, Eq. (13). In Fig. we show the unperturbed rel¬ 
ative velocity found by solving Eqs. (18)-(21| with 
initial relative velocity V^b,o = K-msi and oaryon speed 


of sound Cs = y/iTb/mb- We also plot the same two ve¬ 
locities for an interacting case. All velocities are divided 
by a factor of 1/(2 -I- 1) to eliminate a fiducial redshift 
dependence. 

In the standard case the speed of sound is always below 
the bulk one, which creates supersonic flow of the baryons 
[101 HI]. Including collisions can both raise the thermal 
velocity as well as decrease the relative one, so it reduces 
the Mach number JV to be lower than 1 at lower redshifts, 
which could affect the formation of small-scale structure 


FIG. 9. Values of the relative velocity (solid curves) and the 
thermal speed of sound (dashed curves) divided by (1-1- 2 ). We 
show the collisionless case (black curves) and the case with 
(T 41 = 1 (blue curves), for = 1 GeV. 


function of the dark-matter mass . We show how the 
result asymptotes for very high and very low , and we 
also compare with the CMB-|-Lya analysis in Ref. m, 
shown in dotted-green curve, which is only valid for large 
m^. 



(GeV) 


FIG. 10. Minimum ag (in cm^) detectable as a function of 
nix ill GeV. In black curve we show the results for the case 
with realistic parameters and in blue curve the one with op¬ 
timistic parameters. In dotted-green curve we display the 
current CMB constraint (only valid for ^ GeV). 


VI. CONCLUSIONS 


• Finally, throughout the text we have quoted results 
for rrix = 0.1, 1, and 10 GeV. For lower masses, the 
result is independent of mass, and for higher masses it 
depends on aojirix- We show a larger range of dark- 
matter masses in Fig. |10[ where we plot the minimum 
tJo one could detect at a signal-to-noise ratio of 1 , as a 


We have shown that adding interactions between dark- 
matter and baryons (in a velocity-dependent way a = 
aov~'^) can dramatically change the behavior of the 
baryons during the dark ages. Intuitively it would seem 
that coupling the two fluids will tend to decrease the tem¬ 
perature of the baryons in favor of the dark-matter tern- 
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perature. We have proven, however, that there is an extra 
heating term that appears due to the relative velocity be¬ 
tween the two fluids, which tends to convert initial kinetic 
energy into thermal energy. For a wide range of non-zero 
initial relative velocities and dark-matter masses > 1 
GeV, we find that the heating dominates over the cooling 
of the baryons (if the initial velocity is zero, we find that 
there is only cooling, as expected). 

Heating up the baryons affects the physics of the dark 
ages significantly. A higher baryon temperature will 
translate into a higher spin temperature during the dark 
ages, and hence a smaller brightness temperature of the 
21 -cm line (if one cools down the baryons, the effect is 
the opposite). The dependence of the heating on the ini¬ 
tial relative velocity makes the brightness temperature 
depend on the position in the sky, and hence creates an 
extra source of perturbations to the brightness tempera¬ 
ture, sourced by the relative velocity perturbations. We 
have calculated the power spectrum of these perturba¬ 
tions and compared it to the standard primordial one. 

To find constraints to dark-matter-baryon interactions 
we have studied two probes. First, the global signal dur¬ 


ing the dark ages gets modified by the interactions, and 
assuming an experiment that could detect the global sig¬ 
nal at peak redshift with an accuracy of 1% we can fore¬ 
cast (Jo < 10“^^ cm^. Second, the interactions create a 
new contribution to the power spectrum. We studied the 
case of a realistic ground-based interferometer focusing 
on the end of the dark ages (z = 20 — 30) and found that 
a cross section of (Tq ~ 10 “^^ cm^ could be marginally de¬ 
tected. We also considered a more futuristic experiment 
and found that the minimum cross section one could mea¬ 
sure is (Tq ~ 10 “^^ cm^, more than 2 orders of magnitnde 
better than can be achieved by CMB-kLyo analysis, and 
with a broader mass range. 
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Appendix A: Heating rate 


The baryon heating per unit time is given by 

^ = {r^X^rrib) J J (l^x “ ^b\) Wx “ Vf,| [vcm • (vfc - v^)] , (Al) 

we again resort to writing everything in terms of Vm and w_, this time defining them to be thermal, so they do not 
depend on and we use the general expression a = (Jon", so that 

^ (|V,b +v_|)"« |vc„ ■ (V,b + v_)|, (A2) 

and we can calculate the scalar product term using vcm = + bv- + with a = m^f(m^ + mt) and 

b = (T^ — Tb)/[u^^^{m^ + mb)]. Then the scalar product will be (vcm • v_) = aV^b^+^v^ -|-(a-|-6)V^b'v_ J-v^- (... )■ 
We will have to integrate over Vm and n_, which makes it obvious that the factors proportional to Vm will cancel out, 
as 


d Vmfm — 1 

d * A = 0 . 


(A3) 

(A4) 


Therefore 


so we can calculate the heating rate in terms of two integrals, 

dQb ^bPx^O 


^ = (Z+2 ,) ■ 

ds, 

[aJi(n) J- 6/2 (n)], 


dt (m^ + nib) 


defining 


T/^ + 6 pOO pi 

di (n) = , 3 / / dy {l + + 2xyY'^~''^ {1 + xy), and 

{2Try/^u^^ do J-i 


T/^+6 poo pi 

hin) = 3 / / dy (l-b -b a;(a;-b j/), 

(/TT) / Jo d-l 


(A5) 


(A6) 


(A7) 
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with X = V-jV^b and r = l^-b/uth- We can perform the two different y integrals to find, 

. 2 s(n+i)/ 2 , , (x - l)(n + x + 4) |a: - + (n - a; + 4)(a; + 1)"'+'^ 

/ dy{l + x^ + 2xyy {1 + xy) = ^ ^ 'w ^ ^-, and 

7_i (n + 3)(n + 5)a; 


J dy (^1 + x^ + 2a;?/) {x + y) = 


{x + l)"+‘*[(n + 4)a; — 1] — (a — l)^[(n + 4)a + 1] \x — 1| 
(n + 3)(n + 5)a2 


n+1 


(AS) 


which means that 


h{n) = 


T/n+6 


hin) = - 


T/n+6 

^Xb 


(27r)i/2, 


?J-oo (n + 3)(n + 5) 

r , -.V/2 (^-lf[(^ + 4):r + l]|a-ir^ 

i-oo (n + 3)(n + 5) 


(A9) 


'th a- 


These functions can be expressed in terms of hypergeometric functions, but since we are interested in the n = —4 
case, let us solve just for that 


/i(-4) = - 


/2(-4) = 


V\ 

Xb 


(27r)i/2, 


dxe-^"^"l^x^ 


la- 11 


1 


2e *2 r — v^erf 


'th a-oo 


v\ 

xb 


(2»)i/2t 




(a-1) (27r)i/2Mj|^ 

1 


and 


th -'-oo 


(a-l) (27r)i/2yjj^ 


2e"T. 


We finally find. 


dQb n?hp^(To 


dt (m^ + nib) (27r)i/2ujij 


2{b-a)e +a^erf(^^^ 


or, plugging b and o. 


dQb '^bPx^O 

dt {m^ + mb)'^V^Uth 


,T^-Tb^_4 _ Fir) 


-e 2 _ jYi 


nh 


which in the V- 3> limit (corresponding to r —>■ 0), yields 

dQb mbPxCTo 1 mbPx(To 1 '^x , - T 1 

dt + TOh) (27r)i/2??th + TOb)^ (27r)i/2 nib J ^ ^ 


(AlO) 


(All) 


(A12) 


(A13) 


This expression does not contain a temperature-independent heating (as expected), and matches exactly Ref. [13] for 
n = —4. However, Eq. (A12) works for all r. In the opposite limit, where ^ V- (and then r —> oo), we find that 
the heating is 


TTl-^ TTh^) P'y^ ^0 


dQb 

dt (m^ + nib)'^V^b^ 

which indeed does not cancel even for =Tb, and ends up being independent of temperature. 
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